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We demonstrate precisely what particle physics information can be extracted from a single direct 
detection observation of dark matter while making absolutely no assumptions about the local velocity 
distribution and local density of dark matter. Our central conclusions follow from a very simple 
observation: the velocity distribution of dark matter is positive definite, f(v) > 0. We demonstrate 
the utility of this result in several ways. First, we show a falling deconvoluted recoil spectrum 
(deconvoluted of the nuclear form factor), such as from ordinary elastic scattering, can be "mocked 
up" by any mass of dark matter above a kinematic minimum. As an example, we show that 
dark matter much heavier than previously considered can explain the CoGeNT excess. Specifically, 
m x < mGc can be in just as good agreement as light dark matter, while m x > m,Gc depends on 
understanding the sensitivity of Xenon to dark matter at very low recoil energies, En < 6 keVnr. 
Second, we show that any rise in the deconvoluted recoil spectrum represents distinct particle physics 
information that cannot be faked by an arbitrary f(v). As examples of resulting non-trivial particle 
physics, we show that inelastic dark matter and dark matter with a form factor can both yield such 
a rise. 



I. INTRODUCTION 

What particle physics information can be extracted 
from a detection of events at a dark matter direct de- 
tection experiment? This seemingly innocuous question 
is riddled with subtleties and uncertainties. The main ex- 
perimental uncertainties relate to separating signal from 
background, and translating observed energy in a de- 
tector to nuclear recoil energy. The main astrophysical 
uncertainties are the local properties of dark matter for 
which direct detection is sensitive: the local velocity dis- 
tribution and the local density. 

The average density of dark matter at a galactic ra- 
dius equal to the Sun is reasonably well known based on 
galactic kinematics (for a recent discussion, see PQ). The 
local density could be quite different; simulations that 
try to quantify the likelihood that the local density dif- 
fers from the average density have been done, e.g., [2J. 
The simulations suggest the local dark matter density is 
within one to several orders of magnitude from the canon- 
ical value of 0.3 GeV/cm 3 , with more uncertainty on the 
lower bound than on the upper bound, depending on the 
likelihood one is willing to tolerate. This implies a min- 
imum uncertainty of one to several orders of magnitude 
on the direct detection cross section. (As we will see, 
there remain several additional sources of uncertainty.) 

The velocity distribution is even more subtle. A pri- 
ori not much is known about the velocity distribution, 
although there has been substantial work to estimate its 
form using iV-body simulations [3 8J. These simulations 
are generally excellent tracers of the average properties 
of the velocity distribution for a galaxy that is like the 
Milky Way. They do not, however, have the level of reso- 
lution needed to determine the local velocity distribution. 
Moreover, while most simulations do not incorporate the 



feedback effects of matter, those that do find interest- 
ing "dark disk" structure [5HT2"]. As has been recently 
emphasized in [13] . even isotropic spherically symmetric 
velocity distributions derived from equilibrium distribu- 
tions of dark matter density result in departures from 
the Maxwellian distribution. Local effects on the den- 
sity and velocity distribution by solar and gravitational 
capture have also been considered in [14] . 

Direct astronomical observations of the local stellar 
neighborhood can be inferred to give information on local 
properties of dark matter. The RAVE survey has con- 
strained the local galactic escape speed [TS]. The values 
lie in a fairly large range, from 460 < v esc < 640 km/s, 
depending on the parameterization of the distribution of 
stellar velocities. Observations have also shed light on 
non-equilibrium local stellar motion, such as the pos- 
sibility of "streams" of stars relatively nearby to the 
Sun. Tidal streams can arise from a disrupting satellite 
dwarf galaxy or star cluster. There has been evidence 
of streams in the solar neighborhood for some time [161 — 
I18j . More recent work with newer stellar surveys, in- 
cluding RAVE, SDSS, etc., suggest streams are present 
[H?H2"l"] (although other work does not hnd evidence for 
local streams oriented along the direction perpendicu- 
lar to the galactic disk [22 ). If a tidal stream of stars 
provides a good tracer of dark matter, this or other lo- 
cal structure could significantly affect dark matter direct 
detection measurements [10l [231433] . Generally speak- 
ing, however, prior work has considered streams or other 
non-Maxwellian structure as perturbations on top of a 
(quasi- ) Maxwellian distribution . 

Our conclusion is that while the average density and 
average velocity distribution at the Sun's galactic radius 
do seem to be determined reasonably well from both ob- 
servations and iV-body simulations, the local dark mat- 



ter density and local dark matter velocity distribution 
remain highly uncertain. Breaking with the vast major- 
ity of prior literature on dark matter direct detection, wc 
consider what particle physics properties can be unam- 
biguously extracted while making essentially no assump- 
tions about the density and velocity distribution of dark 
matter. 

We approach this subject making the following simpli- 
fying assumptions: 

1. One direct detection experiment reports several 
events that can be interpreted as dark matter scat- 
tering off nuclei. 

2. Nuclear recoil events are consistent with scatter- 
ing off only one type of nuclei. The scattering 
process off these nuclei is dominated by either a 
spin-independent or spin-dependent process, but 
not both. 

3. No significant time variation of the rate of events 
is found, and no directional detection is observed. 

Since the main purpose of this paper is point out what 
can (or cannot) be extracted from a recoil spectrum in- 
dependent of local density and velocity distribution, we 
focus on what can be obtained from just one experiment 
(or, equivalently, treat several experiments' results inde- 
pendently). Making the second assumption is, for many 
experiments, not an assumption at all when detectors 
are made of a uniform material (e.g. XenonlOO, LUX, 
CDMS, etc.). Even in experiments consisting of mul- 
tiple nuclei with very different masses (e.g. CRESST), 
it would require a remarkable coincidence to have the 
experiment's proportions-by-mass of different nuclei ar- 
ranged such that dark matter were able to scatter off sev- 
eral of them simultaneously with large signal over back- 
ground. For spin-dependent cross sections, typically only 
a few isotopes dominate. Again, it would be a remarkable 
coincidence to have two separate scattering processes (ei- 
ther both spin-dependent or one spin-dependent and one 
spin-independent) contributing to a nuclear recoil signal 
in a single experiment. 

The third point, that we do not consider time varia- 
tion or directional dependence in rates, is a more serious 
omission. Indeed, time variation is itself an intriguing 
signal, especially given the longstanding annual modula- 
tion observed by DAM A [33]. However, we already have 
our work cut out for us with the comparatively mundane 
time-independent signal. Indeed, it is likely that convinc- 
ing evidence of time variation, or directional detection, 
will require years of exposure. Nevertheless, some work 
on the effect of dark matter streams on direct detection 
has already been done in e.g. [23] [25]. We expect to re- 
turn to this pressing issue of the purely particle physics 
implications of time variation and directional detection, 
independent of density or velocity distribution, in future 
work. 

The upshot of these assumptions is that we consider a 
collection of nuclear recoil events and analyze their con- 



sequences. Given the strong bounds on spin-dependent 
scattering from collider searches [35j , we will mostly con- 
centrate on the case of spin-independent scattering. How- 
ever, our formalism and philosophy applies to both types 
of scattering. 

We will find that some characteristics of the recoil spec- 
trum are highly dependent on the astrophysical inputs, 
while others are totally independent of it. We will, of 
course, extract far less information than is usually pre- 
sented in the literature, but the information we do ex- 
tract is far more robust in the sense that it does not rely 
on any properties of the local dark matter halo. How- 
ever, since the only way we can unambiguously probe 
the dark matter velocity distribution and abundance at 
the Earth's position is through direct detection, we argue 
that anything beyond what we claim here ultimately suf- 
fers from dependency on local astrophysical properties of 
dark matter. 



II. KINEMATICS 

We assume dark matter consists of a set of uncolored 
and uncharged particles Xi with mass m i; velocity distri- 
bution fi(vi), and density pi. The dark matter velocities 
are taken to be in the Earth frame. This is a significant 
departure from conventional analyses. It is much more 
typical to take a galactic frame velocity distribution, mo- 
tivated by simulations or simplifying assumptions about 
the galactic distribution, and boost into Earth frame. For 
us, since our velocity distributions f(v) are specified in 
Earth frame, they can be used directly for direct detec- 
tion scattering. However, since time-independent dark 
matter scattering depends only on the magnitude of the 
relative velocity, there is no way to determine the magni- 
tude in the galactic frame outside of a kinematic range. 
This should be kept in mind when viewing velocity dis- 
tributions later in the paper. 

We will mostly consider examples with just one dark 
matter candidate, although the formalism we develop be- 
low will apply to an arbitrary set of WIMPs. For any 
given dark matter particle x, it is useful to first consider 
the kinematics of the most general direct detection recoil 
event. We consider a collision between x with mass m x , 
moving with velocity v in the Earth's frame, and a nu- 
cleus N with mass mjsr, that is stationary. The result of 
this collision, in the most general two-to-two case, will 
be an outgoing dark-sector particle, x'> and an outgoing 
visible-sector particle, N' . We define S x = m' — m x and 
<5jv = tti'n ~ m N, and work in the approximation that 
S x /m x ,SN/mN <C 1, i.e., the mass differences between 
incoming and outgoing particles are small. While gener- 
ally the nucleus is not changed as a result of WIMP scat- 
tering, there are notable exceptions, such as nuclear ex- 
citation [36, 37 and rDM with neutron emission [38ll39j . 

In the center-of-momentum frame, where v is the rela- 
tive velocity of the two incoming particles, the incoming 
momentum is p = fw and the outgoing momentum is p' . 
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Here we use the reduced mass defined with respect to the 
incoming particles, 



m y + m N 



(1) 



The recoil energy of the collision is Er = q 2 /2m' N with 
p' 2 -2pp'cos9 com . (2) 
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q =p 



independent of their mass (and thus, have several candi- 
dates of different masses with similar abundances, using 
for example the WIMP less miracle [13]). 

The event rate of dark matter scattering [44j , differen- 
tial in Er, is determined by 



dR 

dER 



E 



<rvifi(vi(t)) dE ^ , (7) 



The recoil of energy Er, velocity v and cos 9\ a b are related 
by, 
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Define S = 5 X + Sn- If <5 > 0, we can safely perform an 
expansion in 5 jra < 1 to obtain 

1 / m N E R 



^2m N E R \ fx 



(4) 



which taking 8n — > is the well-known result for in- 
elastic dark matter (iDM) [4TM2] . By "safe" we mean 
that our upper bound on w m i n , which is in the far non- 
relativistic regime, automatically implies \S\ <C m x ,m^ 
to allow scattering to be kinematically possible. 

Up to higher order terms in S/m, we obtain an expres- 
sion for the recoil energy 



m N 



E 2 R + 2E R ^(S - [iv 2 cos 2 lab ) + 
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^ -S 2 =Q 



(5) 



The recoil energy is unique for a given fixed scattering 
relative velocity v and nucleus recoil angle #i a b and can 
be solved by the usual quadratic formula, 



E R = \(nv 2 cos 2 6 lah -6) 

TUn 

±{pv 2 cos 2 ^b) 1 / 2 ([iv 2 cos 2 6 Uh - 28) 1/2 



(6) 



This result has the well known feature that the smallest 
recoil energies come from maximizing v 2 cos 2 #i a b, corre- 
sponding physically to head-on collisions at the highest 
velocities available. 



where the sum is over different species of WIMPs, ~ 
Am p is the nucleus mass with m p the proton mass and 
A the atomic number. The recoil energy depends on the 
kinematics of the collision, as described above. Given 
our assumption of no significant time variation in the 
rate, f{vi(t)) — > /(«»), and thus we are effectively ne- 
glecting the Earth's motion around the Sun. This is a 
reasonable approximation so long we are probing veloci- 
ties larger than Earth's velocity in the Sun's frame, i.e., 
w mH > 30 km/s. Typically the maximum speed is taken 
to be » mm = w C arth + w csc , the galactic escape velocity 
boosted into the Earth frame. However, w max is ulti- 
mately determined by the (unknown) details of the dark 
matter velocity distribution in Earth frame. 

Given our assumption of no direction dependent signal, 
we can carry out the angular integral in Eq. Q, reduc- 
ing it to a one dimensional integral where we introduce 
the quantitjj^] fi(v) — J dQf(v). The differential rate 
becomes 
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where we have written 
da, 
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(9) 



in terms of the nuclear form factor F^(Er). There are 
several possible forms for the scattering cross section 
(Ji(v,ER), depending on the interaction, 



<J i0 FI(Er) 

<j i0 (v)FI(Er) 

<7 i0 (v,E R ) 



(10) 



III. EVENT DISTRIBUTIONS 

Our basic assumptions consist of assuming the scat- 
tering process is off only one type of nuclei. We will, 
however, remain general with respect to the possibility of 
multiple WIMPs with different masses, abundances, and 
cross sections. One might think it requires a large coin- 
cidence to have several dark matter particles with cross 
sections large enough to produce events in an experiment. 
However, there are well known counterexamples where it 
can be natural to have the abundance of particles to be 



The different forms for a correspond to functional forms 
of known dark matter scattering that contain velocity 
and/or recoil energy dependence. The first possibility, 
a constant independent of v and Er is the well-known 
isotropic (s-wave) cross section that results at lowest 
order in the non-relativistic expansion from many dark 
matter models. 



1 The velocity distribution is normalized such that f d A vf(v) = 1. 
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The second possibility contains a dark matter form fac- 
tor F Xi (following the standard normalization convention 
F Xi (E R = 0) = 1) and commonly occurs in models of 
composite dark matter |45H47| . Our formalism will han- 
dle the factorizable forms, i.e., the first three of Eq. ( 10 1, 



which incorporates the vast bulk of what has been consid- 
ered in the literature. We will not, however, consider the 
cross sections that contain completely arbitrary nonfac- 
torizable velocity and recoil energy dependence [c.f., the 
most general form written on the fourth line of Eq. ( 10 )]. 



We now turn to the question of what can be inferred 
from a signal in direct detection experiments using Q 
without making any assumptions about fx or the dark 
matter scattering cross section a®. We will however, 
make an assumption about the maximum dark matter 
speed, w max , and we will demonstrate how the derived 
dark matter properties depend on this assumption. 



IV. DECONVOLUTED SCATTERING RATE 



Since the scattering rate (|8]) in any given direct de- 
tection experiment is proportional to the nuclear form 
factor, we first factor it out. This leads to a definition of 
a new quantity, 1Z, that we call the "deconvoluted scat- 



tering rate" - deconvoluted of the nuclear form factor, 
1 dR 



K 



F%{E R ) dE R 



Mim N 



dviVif a (Vi)(Ti(Vi,E R ). (11) 



Some overall factors have been buried into a normaliza- 
tion factor, Mi = NTP Xi /(nfm Xi ). While there are im- 
portant uncertainties in the determination of dark mat- 
ter nuclear form factors from nuclear data [H], this is 
not our concern. Errors on the deconvoluted scattering 
rate ought to take into account nuclear form factor un- 
certainties. 

Next, taking a derivative with respect to E R we find 



dTZ 

dEh, 




duiVifnivi) 



d(Ji(vi,E R ) 



dE R 

fil {Vi Er)) ■ (12) 



For arbitrary 2 — s- 2 kinematics (elastic or inelastic), we 
can replace 



dvi 



dE, 



m N E R - ft dj 

4m N ft 2 E R 



(13) 



This is as far as we can go with a general signal from an 
ensemble of WIMPs with arbitrary cross sections. 

For a single WIMP with a factorizable cross section, 
Eq. ( 11 ) can be used to solve for fi(v) (see also [4"5H52"] h 



fl(Vmin(E R )) 



1 

i%E% - /i 2 5 2 Afa (v min (ER))F 2 (E 




-n 



_l dF 2 (E R ) 

E3(E R ) dE R 



(14) 



This result allows us to gain information on the velocity 
distribution of dark matter evaluated at the minimum 
velocity to scatter for a given recoil energy E R . With 



min ^ 77 



< E" 



we 



scattering data over the range E 
obtain information on the velocity distribution f(v) over 
a range of v. v min {E% in ) < v < u min (^]? ax )- 

For an ensemble of WIMPs, Xh without dark matter 
form factors, the inversion result can be written as 



dn 

dER 



= ^2 w i(v,E R )f il (v) 



where the velocity distributions of the WIMPs are 
"weighted" by the factors 

Wi(v,E R ) = -\ (^J - Jl) Mmoiv) (16) 

For an ensemble of WIMPs with form factors, no simple 
closed form can be written. 



V. /-CONDITION 

There is valuable information that can be extracted 
from Eqs. (14 1 and (15). We know the velocity distribu- 



tion of dark matter must be positive for all v, 
f(v) > , 



(17) 



which we call the "/-condition". Using this condition, 
the right-hand side of Eq. (14) must be positive. Simi- 



(15) 

larly the /-condition also places constraints on the terms 



appearing in Eq. ( 15 ) 



Consider the case of single WIMP with standard elas- 
tic scattering without a dar k m atter form factor, 5 = 
and F 2 (E R ) = 1. From Eq. ( 14) we conclude that the de- 
convoluted scattering rate is always a decreasing function 

Of Er. 

A more striking consequence is reached if a rising de- 
convoluted scattering rate is ever observed. Should there 
be a range of data where the deconvoluted scattering rate 
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satisfies dlZ/dEn > 0, this would signal the presence of 
either an inelastic threshold for scattering or a dark mat- 
ter form factor (or both). This can be seen most easily 
by classifying which terms in Eq. ( 14 ) can be positive or 
negative. 

This general characteristic of the shape of the nuclear 
recoil distribution can thus be used to provide an exper- 
imental signal of non-standard dark matter interactions, 
independent of the local velocity distribution and local 
density of dark matter. This observation is one of main 
conclusions of our paper. 

Eq. ( |14| illustrates how the introduction of non-trivial 
particle physics can allow for structure in the recoil spec- 



trum while still satisfying Eq. (17 1, through either an in- 



elastic cross section or one with a non-trivial dark matter 
form factor. One smoking gun for inelastic dark matter 
is a rise and fall of the (deconvoluted) dark matter scat- 
tering rate. Yet, Eq. ( |l"4"| makes clear that a dark matter 
form factor whose functional form has positive powers of 
En over a range of En within the experimental sensitiv- 
ity also can lead to a rise and fall of the deconvoluted 
scattering rate. Since no direct detection experiment 
can probe F£(E R = 0), the functional form of a dark 
matter form factor is essentially unconstrained]^] Indeed, 
Ref. [S31 |M] demonstrated that specific models of what 
they called "form factor dark matter" can reproduce a 
rising deconvoluted scattering rate, dlZ/dEn > 0, like 
inelastic dark matter. It is interesting to consider simple 
functional forms for the form factor F£(Er), and whether 
they can (or cannot) fake inelastic dark matter given an 
arbitrary velocity distribution. We carry out this study 
in Sec. IVTDl 

The /-condition also allows us to see that there is a 
complete degeneracy between a(v) and fi(v). No one ex- 
periment can separate a velocity-dependent dark matter 
scattering cross section from astrophysical structure in 
the velocity distribution. To the extent that multiple ex- 
periments do not have the same <j(v) for different nuclei, 
it should be possible to break the degeneracy through 
multiple observations. 

It is instructive to consider several possible forms of 
fi (v) and how they map into distributions of energies. As 
a simple starting point, consider a velocity distribution 
that can be approximated by a delta function, 



fi(v) = S(v - v ) 



(18) 



This could correspond physically to a local density of 
dark matter that is contained within a stream moving 
with speed vq in the Earth's frame, or a scattering cross 
section with a resonant feature [35] ■ From Eq. (Ill, 



the deconvoluted spectrum will be a constant provided 
Vmin < an d zero otherwise. In terms of En, the rate 



We are not aware of a way to obtain F^(Er) oc E^ for n < —2 
from a microscopic theory. However, since this does not lead to 
dlZ/dEji > 0, this limitation is not of much concern. 



is constant provided En is below a certain threshold de- 
termined by vq, and zero above it. This result is shown 
graphically for elastic scattering off xenon in the top pan- 
els of Figure [T] 

One can build a discretized version of any fi(v) by 
combining a suitable number of delta functions. Ev- 
ery delta function contributes to the energies below its 
threshold En, and zero above it. In the case of elas- 
tic scattering, the result is a deconvoluted distribution 
which falls off at larger energies. Two other simple cases, 
fi(v) = constant for some range of v and f\{v) falling 
linearly in v are also shown in Fig. [I] 

Several further comments are in order. If we imagine 
a spectrum of several delta functions, the lowest velocity 
contributions may result in recoil energies below an ex- 
periment's detection threshold, and therefore obviously 
cannot be observed. Given that the weights of these 
delta functions are also undetermined, this means that 
the normalization of the velocity distribution cannot be 
experimentally determined. 

Experiments may also be limited in sensitivity at for 
some recoil energy En > £;^ xp - max . Contributions to 
the velocity distribution at very high velocities that con- 
tribute to the recoil spectrum in this range would also 
not be fully probed by an experiment. This is another 
source of uncertainty regarding the normalization of the 
velocity distribution. 

The upshot of this is that there is both astrophysics as 
well as fundamental experimental sensitivity limitations 
in interpreting the overall normalization of a given recoil 
spectrum as arising from a given cross section. This lim- 
itation is fundamental - for an arbitrarily large fraction 
of dark matter moving slowly enough in Earth frame, no 
experiment will be sensitive to such nuclear recoils, and 
thus, no experiment can measure a cross section indepen- 
dently of astrophysics. However the results we find here 
are based on the shape of the recoil spectrum and not the 
normalization, and thus can be determined independent 
of astrophysics. 



VI. CASE STUDIES 

We now consider several examples of what can, and 
cannot, be determined from a nuclear recoil spectrum of 
a single dark matter detection experiment. Given the 
discussion above, in all cases the normalization of fi(v) 
(and thus the cross section) is completely arbitrary. 



A. Elastic Dark Matter 

For the case of elastic dark matter (eDM), defined as 
standard elastic scattering (5 — 0) with no dark mat- 
ter form factor (F^(En) = 1), we can determine a lower 
bound on m x provided one has knowledge of the maxi- 
mum possible velocity, u max . From Eq. |5]), we can solve 
for /i as a function of En and the other parameters. For 



5 



/i(v) 



100 200 300 400 500 600 700 

v [km/s] 



Rate 






100 200 300 400 500 600 700 20 40 60 80 100 

v [km/s] E R [keV] 



FIG. 1. Some examples of the relationship between velocity distribution (LH plots) and observed recoil spectrum, dR/dEn, 
(red dashed in RH plots) and deconvoluted spectrum, 1Z, (blue solid in RH plots). The dark matter mass is taken to be 100 
GeV with elastic scattering off xenon. 



elastic scattering, the result is particularly simple, 



/ m N E R 

^max 



(19) 



demonstrating that the strongest lower bound on the 
dark matter mass comes from the highest recoil energy 
events at the maximum dark matter velocity (in Earth 
frame). We illustrate this bound in Fig. [2j by showing 
the bound on m™ ln as a function of v maK for four possible 
values of the maximum recoil energy from a distribution 
of events where dark matter scatters off xenon. In the 
next section, we will see that the analogous constraints 
on m™ n for inelastic dark matter depends on the inelas- 
tic threshold. 



The deconvoluted scattering rate, Eq. (14), takes on 
the simple form in eDM, 



fx(v) 



AE R dTZ 
m 2 N Afa (v) dE R 



(20) 



The positivity of f(y) (and ao(v)) means that for elastic 
scattering the spectrum of recoil events must be a mono- 
tonically decreasing function of energy. If this is not ob- 
served, we can immediately rule out elastic scattering, 
completely independently of any assumptions about how 
it couples to nuclei or how it is distributed in our galaxy. 

We now discuss what can be determined if indeed a 
falling spectrum is observed. As a surrogate for experi- 
mental data, and to demonstrate our technique, we gen- 
erate pseudo-data for a 100 GeV WIMP elastically scat- 
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FIG. 3. Velocity distributions derived from pseudo-data, with 
5 keV < E R < 80 keV, generated for a 100 GeV WIMP elas- 
tically scattering of xenon. The data is generated assuming a 
Maxwellian distribution for the dark matter, dashed red line. 
For m x = 40, 100, 500 GeV (thick line segments, from bottom 
to top) the derived velocity distributions, fi(v), are shown. 



FIG. 2. The kinematic bound on the dark matter mass, as- 
suming elastic scattering off xenon. Different curves are plot- 
ted for different assumptions about what events are assumed 
to be the highest energy nuclear recoils from dark matter. 
The yellow band corresponds to the typically largest relative 
velocity arising from summing v csc + u cart h < 800 km/s. 



tering off xenon, assuming the standard Maxwellian dis- 
tribution for the dark matter velocity in galactic frame. 
Specifically, our input contains a distribution with char- 
acteristic speed, Wo = 220 km/s (in galactic frame), with 
an escape speed of v esc — 500 km/s (in galactic frame). 
Since we do not consider time-dependent signals, we take 
the Earth's velocity to be a constant, v ear th = 230 km/s 
(in galactic frame). We consider the lower threshold of 
our pseudo-experiment to be recoil energies of 5 keV and 
the upper threshold we take to be 80 keV (the latter 
happens to be below the first zero of the nuclear form 
factor of xenon). These specifications are a reasonable 
approximation of the capabilities of the XenonlOO ex- 
periment [55]. We ignore detector effects such as energy 
resolution and efficiency and assume that these are suffi- 
ciently well known that the recoil energy distribution of 
Eq. can be determined from the data. 



Next, we invert this "data" using Eq. (201 and assume 



the dark matter is scattering elastically. For simplicity, 
we use the same scattering cross section used to generate 
the data. This only affects the normalization of fi(v) 
and does not alter any of our arguments. In Fig. [3] we 
show the derived velocity distributions fi(v) for various 
choices of dark matter mass. As expected, for the correct 
choice of dark matter mass the derived velocity distribu- 
tion agrees with that used to generate the data. However, 
since data is only taken over a finite range of recoil ener- 
gies the velocity distribution is only known over a finite 
range of velocities, corresponding to the v m i n associated 



with the E R , Eq. Q. 

As was discussed in Sec.|ll]it is possible to place a lower 
bound on the dark matter mass by assuming a maximum 
speed for dark matter in our halo, the highest energy re- 
coil events, taken to be 80 keV in the above, then de- 
termine a minimum mass through Eq. Q. However, as 
can be seen from Fig. [3] no such upper bound can be 
made. For all assumed masses used in Fig. [3] the result- 
ing velocity distributions appear a priori to be perfectly 
reasonable: fi(v) is positive and finite. Without further 
mo del- dependent assumptions, or additional experimen- 
tal results (either from another experiment or from rais- 
ing the upper threshold at the first experiment), all dark 
matter masses, m x > m™ m , give a reasonable fit to data. 
Using Eqs. Q and (20) in the limit of very heavy dark 
matter, the derived velocity distribution becomes inde- 
pendent of the dark matter mass and only depends on 
the target. 



B. Elastic Dark Matter: CoGeNT 

A further interesting example of elastic dark matter 
is the recent observation of an excess of low energy re- 
coil events by CoGeNT 56J. The CoGeNT collaboration 
demonstrated that this is consistent with light dark mat- 
ter, 7-10 GeV, where the mass was determined by taking 
the standard Maxwellian velocity distribution. 

Following our procedure above, we have taken a fit to 
the low recoil energy excess at CoGeNT [ST] and reverse- 
engineered the distribution to determine the needed ve- 
locity distribution for (much) larger dark matter masses 
in Fig. [4] It is easy to understand the shift in the veloc- 
ity distributions. For light dark matter, m x < m,Q ei the 
maximum recoil energy is _E™ ax ~ m v t, max/ m Gc- Once 
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FIG. 4. Velocity distributions derived from the CoGeNT 
data. The dashed red line is the ordinary Maxwellian ve- 
locity distribution for an elastically scattering dark matter 
particle that fits the CoGeNT excess with a mass of 8 GeV. 
The derived velocity distributions, fi(v), are shown for m x = 
8, 20, 100, 500 GeV (thick line segments, from bottom to top). 
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FIG. 5. Predicted recoil spectra at a Xenon experiment for 
the velocity distributions that fit the CoGeNT excess shown 
in Fig. [4] The lines correspond to m x — 8, 20, 100, 500 from 
left to right. Notice that the normalization corresponds to 
thousands of events at XenonlO or XenonlOO. 



the dark matter mass is larger than about mce, the maxi- 
mum recoil energy asymptotes to Ef{ ax ~ mG C ^ 11H . The 
shift in w max from m x ~ 7 GeV to m x ^> tog is thus 
a factor of ~ 10. Indeed, from Fig. [4] we see that the 
largest minimum velocity shifts from about 500 km/s to 
about 60 km/s as the mass is increased from 8 GeV to 
500 GeV. 

By itself, the CoGeNT data can thus be fit by any dark 
matter candidate above the kinematic minimum, which 
for v max ~ 800 km/s is about 6 GeV. The conclusion 
that CoGeNT implies "light dark matter" is thus not 
warranted if the local dark matter density and velocity 
distributions can be freely adjusted. 

The obvious objection to the large dark matter, small 
velocity distribution dark matter interpretation of the 
CoGeNT data is that other direct detection experiments 
should be sensitive to larger masses, and potentially rule 
it out. We can test this assertion by translating the Co- 
GeNT observed energies into nuclear recoil energies (fol- 
lowing [58]). and then compute the predicted spectrum 
at an experiment using xenon. This assumes, of course, 
that whatever scattering process is occurring in germa- 
nium also occurs in xenon. 

In Fig. [5] we show the predicted recoil spectrum at 
a xenon experiment, for each of the velocity distribu- 
tions show in Fig. |4j In each case, the prediction is large 
numbers of recoil events below 6 keVnr. This may or 
may not be allowed by existing data from XenonlO [ST?] 
and XenonlOO [55) . The principle difficulty is determin- 
ing the sensitivity of these experiments to very low re- 
coil energies, specifically the conversion factor L c r{Er). 
There has been a spirited discussion on this point [60F 
[66] . Clearly, masses much heavier than considered by 
CoGeNT are allowed, while arbitrarily large masses (rep- 
resented by 500 GeV) may be constrained by this existing 



or future xenon data. In any case, our central conclusion 
is that qualitatively larger dark matter masses can fit the 
CoGeNT data if the velocity distribution is adjusted as 
we illustrated above. 



C. Inelastic Dark Matter 

Due to the splitting of inelastic dark matter, v m i n is 
not monotonic in Er, which means that the highest en- 
ergy recoil events may not give the strongest bound on 
the dark matter parameters. The expression for v m i n , 
Eq. Q, has a minimum, and correspondingly the decon- 
voluted recoil spectrum, 1Z, has a peak at 

£P eak = S , (21) 

m x + tojv 

the observed spectrum, dR/dEji, typically has a peak at 
a lower energy due to the nuclear form factor. Whether 
the highest or lowest energy bins of an experiment place 
the strongest constraints on m x and S depends on where 
this peak falls. 

If i?^ cak is large enough there may even be an observ- 
able gap between an experiment's lower threshold and the 
first recoil events - this bump-like spectrum aids in fit- 
ting the energy spectrum of the DAMA modulation data. 
On the other hand, if E R cak is below the lower threshold 
it is not possible to distinguish iDM from conventional 
elastic dark matter. In Fig. [6] we illustrate these points 
for the case of a xenon experiment which records data 
over the range 5 keV < Er < 80 keV. 

Unlike eDM, where the recoil spectrum must be mono- 
tonically decreasing, the mass splitting of iDM allows for 
the spectrum to be a rising function for Er < E R ciik . 
This cannot be mimicked by elastic scattering without 
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FIG. 6. The astrophysics independent allowed region of iDM parameter space (unshaded region) for a recoil spectrum with 
events, in a xenon detector assuming a maximal speed for dark matter of 500 km/s (LH plot) and 800 km/s (RH plot) in the 
Earth's frame. In both plots the two curves assume events observed with 5 and 80 keV. The lower (red shaded) region is 



model independent and is the region in which the peak, Eq. (21 1, in the deconvoluted recoil spectrum lies below 5 keV and 



iDM cannot be distinguished from eDM. With reliable knowledge of fi(v) this region extends further upwards. 



a dark matter form factor, no matter what (physically 
allowed) form the velocity distributi on takes. Further- 



more, at i??, cak the denominator in (|14|) has a zero but 



since fi(v) must be finite the recoil spectrum must have 
a maximum at this same energy so the numerator goes 
to zero as well. Thus the observation of a maximum in 
the recoil spectrum determines £^ cak and reduces the 3 
dimensional iDM parameter space by one dimension, re- 
lating m x and S, independent of knowledge of the dark 
matter velocity distribution. 



D. Form Factor versus Inelastic 

Although the spectrum of iDM cannot be faked by sim- 
ple elastically scattering DM, regardless of what physical 
velocity distribution is used, it may be faked by elasti- 
cally scattering dark matter that has a non-trivial form 
factor. Due to the unique kinematics of iDM, induced 
by the splitting of the states, a deconvoluted iDM spec- 
trum has a gap in the low energy spectrum with no recoil 
events. The size of this gap, and thus the position of the 
first events is 



pgap _ 



m N 



v — 



(22) 

This can be obtained directly from Eq. |7| by calcu- 
lating the lowest recoil energy for the highest velocity 
(v = v max ) for a head-on collision (cos#i a b = 1). The 
spectrum then increases to £^ oak , Eq. (21), before again 



decreasing. Near the peak, the spectrum is well fit by a 
power law but away from the peak the spectrum becomes 
more exponential in nature. This increase in 1Z between 



Ef^ p and £^, cak would require an unphysical fi(v) < if 
the dark matter has no form factor. However, it is pos- 
sible that a non-trivial F x (Er) could "overpower" this 
increase and allow iDM spectra to be fit by eDM with a 
form factor, albeit for a non-Maxwellian velocity distribu- 
tion. We demonstrate this below for the simple example 
of a form factor that is a single power in En, F x ~ E^, 
for a sufficiently high n the resulting f\{v) will be sensi- 
ble. 

Although the increase in the recoil spectrum can be ac- 
commodated by elastically scattering dark matter with a 
non-trivial form factor, the existence of a gap with iDM 
cannot be faked with a simple power-law form factor. 
This is because the threshold to up-scatter results in 
a step function turn on in the spectrum, and no such 
step function results from a simple polynomial form for 
F*(E R ). 

To illustrate the ability of iDM to be faked by form 
factor dark matter. We consider the simple case of F x = 

\<t) w ith the normalization qo — 10 -3 /i and we take 
n G [—1, 10] and investigate whether there is a physically 
sensible velocity distribution that allows form factor dark 
matter (with the same WIMP mass) to fake iDM with 
m x = 100 GeV, with a Maxwellian velocity distribution, 
for various splittings 5. This region in which iDM can 
indeed be faked is shown as the shaded region in Figure [7] 
and an example velocity distribution that achieves this 
is shown in Fig. [8] 



VII. DISCUSSION 

We have shown that certain particle physics proper- 
ties can be determined from the nuclear recoil spectrum 
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FIG. 7. In the shaded region elastically scattering dark matter 
with a form factor F x oc q 2n can fake iDM with splitting S. 
The dark matter mass was taken to be 100 GeV, scattering 
off xenon, and the iDM spectrum was assumed to come from 
a Maxwellian distribution. 



of a dark matter direct detection experiment independent 
of the local astrophysical density and velocity distribu- 
tion. These properties are most easily uncovered from 
the deconvoluted scattering rate, that we call 7Z. We 
advocate that experiments present their recoil spectrum 
data in terms of the deconvoluted spectrum, from which 
the shape of 1Z as a function of En can be read off inde- 
pendently of the nuclear form factor. 

The kinematics of scattering leads to a lower bound 
on dark matter mass, given a maximum speed for the 
WIMP in the Earth's frame, independent of the details 
of the dark matter velocity distribution. We have also 
shown that an upper bound on the mass of dark mat- 
ter, however, cannot be determined independently of the 
density and velocity distributions of dark matter. This 
was demonstrated for both "fake" Xenon data, but also 
for the CoGeNT excess. Perhaps the most striking con- 
sequence is that we find the CoGeNT excess can be fit 
by dark matter with masses much heavier than has been 
generally considered. Much of this parameter space, up 
to m x ~ togc, is no more strongly constrained than light 
dark matter, given an appropriate choice of velocity dis- 
tribution. For m x > rriQ e , the CoGeNT data may or may 
not be allowed by existing XenonlO and XenonlOO data, 
which ultimately depend on the details of L c g(Eji) and 

Qy 

Our analysis comparing different masses, as well as 
comparing and contrasting inelastic dark matter with 
form factor dark matter, was done without any exper- 
imental errors. It would be interesting to carry out sim- 
ilar analyses with experimental errors. For example, the 
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FIG. 8. Top: a recoil spectrum, at a xenon based detector, 
for iDM with m x = 100 GeV and S = 35 keV within the 
range 5 < E R /keV < 80. Bottom: the iDM Maxwellian 
velocity distribution is shown as the the red dashed curve. 
The solid gray curve corresponds to the velocity distribution 
that exactly matches the recoil spectrum shown in the top 
figure, but with dark matter, of the same mass, scattering 
elastically with a form factor F x oc q 4 . 



so-called "gap" in inelastic dark matter may not be re- 
solvable due to experimental errors, so that it may be 
reproduced by a smooth dark matter form factor. 

Breaking from the vast majority of past literature, we 
do not impose any requirements of the velocity distri- 
bution of dark matter, except for assuming there is an 
upper bound (that might or might not be equated to 
Woarth + v csc in Earth's frame). In fact, our velocity dis- 
tributions were taken to be in Earth's frame, and thus 
to compare with most other literature one would have to 
boost into galactic frame. 

Interestingly, however, there are potential physical re- 
alizations of some of the more unusual velocity distri- 
butions found in this paper. In particular, distributions 
with velocities that are smaller than the Earth's velocity 
in the galactic frame (~ 230 km/s) could be interpreted 
in at least two ways. One is dark matter that streams 
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in a direction parallel to the motion of the Sun around 
the galaxy. This could be a local stream, or could be 
a "dark disk". Another interpretation is that the lower 
velocities arise from dark matter that is weakly bound 
to the Earth/Sun/ Jupiter system. In both cases, relative 
velocities in the tens, rather than hundreds, of km/s are 
physically realizable. 

There are still many questions and many inter- 
esting avenues for future exploration of astrophysics- 
independent implications of dark matter. Time- 
dependence in the scattering rate is perhaps the most 
important, given the ongoing interest in the annual mod- 
ulation observed by DAMA. In the future, directional 
dependence will also provide an excellent way to deter- 
mine the nature of the dark matter velocity distribution 
we find ourselves in. Understanding the particle physics 
implications of these observations is left to future work. 
With information about dark matter from other sources, 
such as mass determinations from the LHC, the approach 
presented here could allow the interpretation of direct 
detection results as a direct measurement of the local as- 



trophysical properties of dark matter. Finally, we have 
not considered the implications of multiple nuclear recoil 
spectra from different experiments. This too is an very 
interesting and timely subject that we hope to see work 
on very soon [ST] . 
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